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_i , It is shown that the famihes of generaUzed matrix ensembles recently considered 

O ' which give rise to an orthogonal invariant stable Levy ensemble can be generated 



by the simple procedure of dividing Gaussian matrices by a random variable. The 
nonergodicity of this kind of disordered ensembles is investigated. It is shown that 

"j^ , the same procedure applied to random graphs gives rise to a family that interpolates 

c/3 , between the Erdos-Renyi and the scale free models. 



^ i The classes of random matrix ensembles introduced by Wigner in the 50s have found 

' ' a great sucess partly after being connected with quantum manifestations of chaos in 

^ ■ physical systems [ll]. In turn this success generated a great activity and extensions and 

P ! generalizations of those ensembles have occurred. In obtaining the Gaussian ensembles, 

Wigner adapted the Wishart ensembles well known to statistitians. Some of the exten- 
sions of the Gaussian ensembles can also be considered as applications of known processes 
<^ • in statistics. For instance, models to describe symmetry breaking have been constructed 

G^ ! by adding two random matrices, one block diagonal and the other its complement [2] . 

Here we consider a random process in which a new random quantity is generated by 
cn ' taking not the sum but the ratio or the product of two other independent ones. 

In a previous paper [3(], an alternative to Shannon information entropy, namely Tsallis- 
Renyi information [^ was used to introduce a new family of generalized matrix ensembles 
(see also |5|). One of the main features of this ensemble is the power-law characteristic 
of its statistical properties. In particular, it was shown that individual matrix elements 
- behave like the elements of the so-called Levy matricesfO"] (after the publication of Ref. 

1^ ! [Sj, Klauder and Muttalib obtained an even more general family [7] on similar lines), 

cd I One of the purposes of this note is to show that all these families can be obtained, in 

fact, by the following simple procedure. Let Hda) be a random matrix of dimension A^ 
and variance l/2a^ and let its probability distribution be 

//2 



Pg{H; a)= i — \ exp {-aptiH'^) 



The matrices of the Gaussian ensemble are specified by a. In ([T]), / is the number of inde- 
pendent matrix elements f = N -\- (3N{N — l)/2 and f3 is the Dyson index /3 = 1, 2, 4 for 
GOE, GUE and GSE (here and in what follows the subindex G indicates Gaussian). The 
distribution is normalized with respect to the measure dH = Yli dHu Y[j>i Y[k=i V^dHf-. 
Take now a positive random variable ^ with a normalized density probability distri- 
bution w{^) with average ^ and variance o"| and introduce a new matrix ensemble by 
the following relation (product of random variables has been considered in the context 



of covariace matrices^ 

In this way, an external source of randomness is superimposed to the fluctuations of the 
Gaussian matrix Hcia)- A random process in which there is a competition between two 
types of random variables is typical of disordered systems or, in the case of Ising models, 
spin glasses|9|. As the two types of randomness are independent, one can be kept frozen, 
quenched in technical terms, while the fluctuations of the other continue to operate. 
Here the disorder is represented by ^ which is the quenched variable in opposition to the 
randomness of the Gaussian matrices. We may refer to ([2]) as a disordered ensemble. 

From ([2]), we deduce that the joint distribution of a set of n < / matrix elements is 
given by 

where hi = Hij for the diagonal and hi = V^Hij for the off-diagonal elements. Eq. ([3]) 
shows that matrix elements are correlated. As a particular case, for n = f, (^ leads to 
the ensemble distribution 

PiH- a)= I d^wiO (~) '" exp f-^tr//^") (4) 



< ) \ i 

where the term after w(^) is just ([T]) with a replaced by ot^/^- Expressions like (jl]) are 
being considered as instances of superstatistics[10]. 

The relation ([2]) makes straightforward to do numerical simulations in terms of Gaus- 
sian matrices. However, it may also be useful to directly generate matrices of the ensemble 
(taking into account the corrrelations among their elements). This can be done through 
the identity 

p{h,,...,hf)=p{h,)\{^^hl^ (5) 

where each fraction gives the conditional probability for the rath element once the n — 1 
previous ones are given. This equation provides a way to sequentially generate all the 
matrix elements. At each step, a new element, say the nth, is sorted using Eq. ([3]) that 
implies 

where he is a Gaussian variable and C,n is another random variable sorted from the 
distribution 




wnio = ^(oe"-'^/'exp -V y M]/ rfe^(oe^"-^^/'exp 




which is univariate since all the previous n — 1 elements have already been determined. 

By generating matrices fixing, in the process, a set of values ^1,^2, •••)^/ we are, in 
the language of the disordered systems, quenching the disorder. The differences among 
matrices generated with different sets of ^ depend on the width of the distribution w{C,) 
and one can expect that for wide w{^) the large spread among the matrices will give rise 
to a nonergodic behavior. 

Turning now to eigenvalues and eigenvectors, we observe that we have an ensemble 
invariant under unitary transformation in which, as it occurs with the Gaussian ensem- 
bles, the joint distribution of eigenvalues and eigenvector factorizes. The eigenvectors 
behave as those of the Gaussian ensembles and we can integrate them out to obtain for 
the eigenvalues the joint distribution 



where Xi = \/a^/^Ei and 

N 
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Pg{xi,...xn;-) = K^^exp -^z^xl j JJ|xj-x/, (9) 



fc=l / j>i 



with K]^ being a normalization constant. 

From ([H]), measures of the generalized family can be calculated by weighting the cor- 
responding measures of the Gaussian ensembles with the w{^) distribution. Integrating 
for instance ([8]) over all eigenvalues but one and multiplying by N, the eigenvalue density 



is expressed in terms of the Wigner's semi-circle law|ll 



as 



p{E;a) = ^ j diw{0{i/if^^f^N~h^^iE^l (10) 

where the condition a^E"^ < N on C, has to be satisfied. 

As previously stated, the introduction of the disorder represented by the vari- 
able ^, breaks in principle the ergodicity of the Gaussian ensembles. Let N{L) = 
Ie~l 2 dE'p{E') be the average number of eigenvalues in the interval [E — L/2, E + L/2] 
for an ensemble with eigenvalue density p{E). The variance S^(-£/) of the number of eigen- 
values in that interval can be expressed in terms of the two-point correlation function 
R{E^,E2) hy 

l'E+L/2 i'E+L/2 

T?{L)= dE, dE2R{E^,E2) + N{L)~N\L). (11) 

Je-L/2 JE-L/2 

Ergodicity implies 1 121] the vanishing of 

Varp = [p{E)fT?{L)/L^. (12) 

when L —>■ 00. For the disordered ensemble we have 

S^ (L) = I diw (0 [S^(L) - Ng{L) + NliL)] + N{L) - N\L). (13) 

with Ng{L) calculated with the Gaussian density. In (fTSll . nonergodicity will result if 
the quadratic terms do not cancel. Indeed, in this case, a parabolic contribution for 



large L survives and the variance of the density fluctuations given by Eq. flT^ does not 
asymptotically vanish. 

Consider now a particular choice of the distribution w(^). Note that the factor multi- 
plying the Gaussian matrices in Eq. ([2]) acts on the variance of the Gaussian ensembles. 
In order to investigate ensembles showing heavy-tailed densities it is convenient to choose 
w{^) to be the gamma distribution 

w{0 = exp(-oe^""Vr(0 (14) 

that becomes a x^ distribution for integer 2^. From (IT4l) a^ = y f , showing that ^ controls 
the behavior of the distribution w{C,). It becomes more localized when ^ increases and we 
should then expect to recover the Gaussian ensembles. However, for smaller values of ^, 
departures from the Gaussian case will be observed. Indeed, by substituting (TT^ in (jlj) 
we find 



/ 



r A 



Pa y [—ij ( j3a^ ^\^-o 

<V r(e) V^ e 



P{H;a,0 = {^] -^}^{l + ^tTH'] (15) 



for the ensemble density distribution, where 



1 - f 

^ + ^, with g> 1. (16) 



q-1 " 2 

Eq. (IT^ is just Eq. (4) of |3|. In [3|] it was derived using a generalized maximum entropy 
principle [4l] with q being identified with the Tsallis entropic parameter. 



Substituting ([11]) in (|3]) for n = 1[13 



,, ,, //?«V r(e + i/2) f^ ^ Pa,,Y'^-"' ,,^, 

p(/^;«,0=^^j r(^-) i^ + T'J ^''^ 

for the density distribution of a given matrix element. Since for large \h\ , ppih] a,^) ~ 
1/ 1^1 ) (I17p exhibits the power-law character of the distribution. It is important to 
remark that, apart from the lack of independence, the marginal distribution of the matrix 
elements have the same kind of distribution, namely one with an asymptotic power-law 
behavior, as the i.i.d. ones of the ensemble of Levy matrices[6j 

In Fig. 1 the eigenvalue density for three realizations of the ensemble generated 

using the above random process with ^ = 1/2 is histogrammed and compared with the 

1 1 



l+x^ 



semi-circle law. We recall that for ^ = 1/2 the matrix elements are Cauchy, - 
distributed(see Eq. flT7|) ). It is seen that the individual matrices of large sizes are 
Gaussian ensemble matrices as they should. As a comparison, in Fig. 2, it is shown 
the eigenvalue density of just one Levy matrix of large size whose matrix elements also 
follow the Cauchy distribution. We can see that although individual matrix elements 
of the two ensembles are identically distributed, their eigenvalue density behaves in a 
completely different way. While individual Levy matrices of large sizes do not depart 
from the ensemble average, matrices generated according to ([6]) show large fluctuations. 
Of course, the result shown in Fig. 1 indicates strong nonergodicity. This is confirmed 
by the ensemble number variances shown in Fig. 3. The parabolic behavior seems to 



persist even for large values of the parameter ^, showing that the ensemble is noner- 
godic. Consequently, averages performed running along one spectrum do not coincide 
with averages over the ensemble of matrices. 

Other systems in which nonergodicity may play an important role are networks and 
their associated graphs. We now show how the present approach can be applied in 
random graph theory |l4|. A graph is an array of points (nodes) connected by edges. It is 
completely defined by its adjacency matrix A whose elements Aij have value 1(0) if the 
pair (ij) of nodes is connected (disconnected). The diagonal elements are taken equal to 
zero, i.e. An = 0. Adjacency matrices of graphs in which the connections are randomly 
set, are real symmetric random matrices. The classical random graph model proposed 
by Erdos-Renyi (ER) is simply defined by giving a fixed probability p that a given pair 
of nodes is connected, independently of the others |l5l|. 

We start by showing that the ER model can be considered as the equivalent in random 
graph theory to the Wigner model of Gaussian matrices. In fact, the joint matrix element 
distribution of its adjacency matrix A can be written as 

PERiA,a) = [1 + exp(-a)]-^exp (-ftrA^) (18) 

where / = ~ with N^ the size of matrix, being equal to the number of nodes. Eq. 
( fTSj) is just the defining equation ([T]) of the GOE (/3 = 1) ensemble with the constraint 
that the matrix elements can only take the values and 1 imposed by the measure 

N 

dH = \[ dHu5{Hu) n ^dH,j [6{H,j) + 6{1 - Hi,)] . (19) 

1 j>i 

From flTSl) it follows that the marginal distribution of a given matrix element, say Aij, is 

PEuiA,, a) ^ ^^Rt^ . I T^' 11 ^^^ - I (20) 

1 + exp(-a) [ i+exp(-a) ' if ^ii = 0' 

which means that the probability p that defines the ER model is connected to the pa- 
rameter a by the relation 

a = ln(- - 1). (21) 

p 

Since the probability p is defined in the interval [0, 1], the domain of variation of a is 
]cx3, — cx)[. This suggests that the statistical properties of the ER model must show a 
symmetry with respect to the point a = (or p = 1/2). 

It is important to remark that although Eq. flTSj) has the same structure as Eq. 
([1]) there are striking differences between the two models. Despite the presence of the 
trace in flTSl) . the discrete nature of matrix elements imposed by the measure, Eq. flT9|) . 
destroys the rotational invariance and prevents the factorization of the joint distribution 
of eigenvalues and eigenvectors. The parameter a is just a scaling parameter in the 
Gaussian case. In contrast, the properties of ER model depend strongly on the value 
of the probability p, and here a plays an essential role. Notice also that, contrarily to 
the Gaussian cases, the adjacency matrices form an ensemble with a finite number of 
matrices. It is convenient in the study of the graphs, to introduce the scaling p ~ N~^ 
{z > 0). For instance, connectivity properties of the graph are characterized by z. 



An analytical expression of the spectral density for arbitrary values of the probability 
p and matrix size A^ is an unsolved problem [16]. However, when p is fixed and N is very 
large, the density can be deduced in the following way. A is a symmetric non- negative 
matrix with maximum principal eigenvalue, Ei, its value is close to the nonzero eigenvalue 
of the constant matrix < A> with elements equal to the average of the A-elements, i.e. 
< A >ij= p. As the only nonzero eigenvalue of a constant matrix is equal to the product 
of its size by the element, we conclude that Ei = pN. Because of this linear dependence 
with A^, for fixed p the largest eigenvalue grows faster than the others as the matrix size 
increases. In this case, for very large matrices the other eigenvalues have asymptotically 
the same eigenvalue density of the eigenvalues of the matrix A— < A > . This density 
can be obtained from the moments of the trace of the powers of the matrix and one finds 



that it obeys the Wigner semi-circle law|14l 



f 1 ^ \/ANa'^ - E"^ if \E\ < \/ANa^ 
PERiE,a) = I 2.a2V47va is , II is < ^/Zil^ (22) 

'^^^ I 0, if \E\ > VaN^ ^ ' 

where a^ is the variance of the matrix elements given by 



a'=pil-p) = ^—^^—^. (23) 



1 

4 cosl?(a/2) ' 

The above arg ument fails if p ~ 1/A^ (2; ~ 1) in which case deviations from the semi-circle 
appear [la, llTI . 



P{A;a) = I d^wiO ,^ ':' t ZTT/- (24) 



We now introduce a disordered model of random graphs by defining an adjacency 
matrix with a distribution 

exp (-f trA2) 

[l + exp(-«0]^ 

Therefore this generalized model is a superposition of Erdos-Renyi random graphs with 
distribution P{A,a^) weighted with w{C,) exactly as in (jlj) for the disordered Gaussian 
ensembles. Again the width of the distribution of if (^) is a controling parameter and as 
remarked before the parameter a also plays an essential role. In particular, for a = the 
ensemble is just the ER with p = 1/2. 

From Eq. fl2^ we can derive the probability distribution for a set of matrix elements 
and use Eq. (^ to define a random process entirely equivalent to the one used to generate 
matrices of the disordered Gaussian ensemble. As before, a set of probabilities pn with 
n = 1,2,3...,/ is sequentially generated and, from them, each new matrix element is 
obtained taking into account those already determined. This means that Eq. (12^ defines 
a model of a disordered correlated graph in which new attachments depend on the ones 
already existing. 

As in the case of the Gaussian ensembles, statistics of the averaged graph (our model) 
are averages over the ER statistics. For instance, the eigenvalue density is 



where 



p{E; (^) = IJ^ d^HO cosh{^)^N - cosh'{^)E^ (25) 



cosh-^(^). (26) 



We now make for w{C,) the same choice as before, namely Eq. (1141] . As before we 
expect for large values of ^ small fluctuations around ER, whereas for small values they 
will become large and will govern the asymptotics. 

In Fig. 4 we display the density of eigenvalues of the adjacency matrices. When going 
from z close to 1 to 2; close to 0, the density goes from a highly picked density with 
heavy tails towards a Wigner semi-circle, showing a crossover which is reminiscent from 
a scale-free to an ER graph. 

In summary, we have discussed a new method to introduce matrix ensembles which 
preserve unitary invariance presenting distribution with heavy tails. The price to pay to 
preserve unitary invariance is i) to abandon the statistical independence of the matrix 
elements ii) to abandon the ergodic property (equivalence of spectral and ensemble aver- 
ages). There are cases, however, in which only ensemble averages make sense. Consider, 
for instance, the behavior of individual eigenvalues. Recently, extreme eigenvalues have 
been a matter of great interest due to the discovery that the distributions they follow, 
the so-called Tracy-Widom|18] in the case of the Gaussian ensembles, show universality 
and have wide applications [19]. The same authors have found growing systems in which 
an external source induces the extreme values to have a behavior in which there is a 



competition between their distribution and a Gaussian[20|. In a paper in preparation, 
we show that the disordered ensemble can be a useful model for this kind of systems. 

Let us finally mention that the method discussed here (Eq. ([2]) with the choice Eq. 
(TT4|) for the probability density function w{C,)) was intended to rederive and to give new 
insight on models previously studied. By making other choices for w{^) new models 
preserving orthogonal invariance may be introduced (see also [7|). 

We thank L. Pastur and W. F. Wreszinski for fruitful discussions. This work is sup- 
ported in part by the Brazilian agencies CNPq and FAPESP. 
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Figure Captions 

Fig. 1 The eigenvalue density of three matrices of size A^ = 300 generated using Eqs. 
([6]) and (fT4l) with ^ = 1/2 compared with Wigner's semi-circle law. 

Fig. 2 The eigenvalue density of one Levy matrix of size A^ = 600 whose elements are 
Cauchy distributed compared to a Cauchy distribution. 

Fig. 3 Full lines: the number variances calculated with Eq. (I13p for the values 
^ = 5, 10, 20, 50 and 200 as indicated in the figure; dashed lines: the linear Poisson 
number variance and the GOE number variance. 

Fig. 4 The eigenvalue density of the disordered random graph model calculated with 
Eqs. fl25l) and flT^ with ^ = 1/2 and for values 0.2, 0.3 and 0.8 of the scaling parameter 
z. 
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